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We derive the expression for the nonperturbative coupling of a hadron to two external currents 
from its underlying structure in QCD. Microscopically, the action of each current is resolved to a 
coupling with dressed quarks. The Lorentz structure of the currents is arbitrary and thereby allows 
to describe the hadron's interaction with photons as well as mesons in the same framework. We 
analyze the ingredients of the resulting four-point functions and explore their potential to describe a 
variety of processes relevant in experiments. Possible applications include Compton scattering, the 
study of two-photon effects in hadron form factors, pion photo- and electroproduction on a nucleon, 
nucleon-pion or pion-pion scattering. 
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I. INTRODUCTION 

Pion- and photon-induced reactions with the nucleon 
have been the main source of experimental information 
on the structure and properties of hght hadrons. In addi- 
tion to the data collected in Nt: scattering, pion photo- 
and electroproduction experiments have provided new in- 
sight in the electromagnetic properties of nucleon reso- 
nances such as the A(1232) or A/'(1440) [1-5]. More- 
over, virtual Compton scattering in the low-energy region 
determines the nucleon's generalized polarizibilities and, 
at higher energies, allows to resolve its spatial and spin 
structure [6-8], thereby providing a connection to the 
underlying substructure in Quantum Chromodynamics 
(QCD). The Compton scattering amplitude also deter- 
mines the magnitude of two-photon effects in the extrac- 
tion of nucleon form factors [9]. 

Depending on the studied momentum range, different 
theoretical approaches are employed to describe these re- 
actions. In deeply virtual Compton scattering (DVCS), 
the factorization property of QCD allows to disentan- 
gle perturbative and nonperturbative contributions and 
model the soft parts of the process by generalized par- 
ton distributions [10-12]. Pion- nucleon and pion pho- 
toproduction on the other hand are analyzed by uni- 
tarized chiral perturbation theory [13], phenomenologi- 
cal partial-wave analyses [14-16] or dynamical coupled- 
channel models [17-21]. Here, the dominant s— and 
t— channel poles are usually implemented via potentials 
and the scattering matrices are generated by iterating 
Bethe-Salpeter equations at a purely hadronic level. Fea- 
tures such as unitarity and analyticity of the S*— matrix, 
electromagnetic gauge invariance and crossing symme- 
try are desired but not always realized to full extent, see 
Ref. [21] for a discussion of this issue. Nevertheless, these 
approaches provide important theoretical input to the 
experimental data analyses and are necessary to recon- 
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struct baryon form factors from the experimental cross 
sections. 

A desirable alternative would be a nonperturbative cal- 
culation of such scattering amplitudes within QCD it- 
self, i.e., by resolving the underlying dynamics in terms 
of quarks and gluons. Typical diagrams that a micro- 
scopic description should recover in various kinematical 
ranges are shown in Fig. 1 for the case of pion photo- 
and electroproduction. Microscopically, the pion and the 
photon couple to dressed quarks. This leads to handbag 
diagrams (Fig. la) which, in the analogous case of vir- 
tual Compton scattering, become dominant at large pho- 
ton virtualities, as well as so-called cat's-ears diagrams 
(Fig. lb) that comprise two active quarks. However, 
the relevant contributions in the low-energy regions are 
others, namely processes that involve s— channel nucleon 
resonances and i— channel meson exchange (Fig. Ic and 
Id). A comprehensive description should resolve those 
diagrams in terms of their quark and gluon substructure 
as well. The purpose of this paper is to derive such ex- 
pressions in a systematic way. 

We work in the framework of Dyson-Schwinger equa- 
tions (DSEs) of QCD [22-24]. They interrelate QCD's 
Green functions and thereby provide access to nonper- 
turbative phenomena such as dynamical chiral symmetry 
breaking and confinement. In combination with covari- 
ant bound-state equations, i.e., Bethe-Salpeter equations 
(BSEs) and Faddeev equations, they present a compre- 
hensive framework for studying hadron properties from 
QCD. The approach allows to probe the substructure of 
hadrons at all momentum scales and all quark masses. It 
has been applied to compute a variety of hadron observ- 
ables such as meson and baryon masses, wave functions, 
leptonic decay constants, elastic and transition form fac- 
tors; see [24-28] and references therein. 

Building upon that, a systematic and nonperturbative 
procedure to construct hadron matrix elements from the 
underlying Green functions is the 'ganging-of-equations' 
method [29-31]. The basic idea is to couple an external 
current such as a photon to all internal building blocks 
to arrive at a quark-level description of hadronic matrix 



Vi" 






(a) 



(b) 



N,N*,A,... / \ n,p,. 

(c) (d) 



FIG. 1. (Color online) Various contributions to the four-point function, here in the case of pion photoproduction on the nucleon, 
that must be recovered in a microscopic description: handbag (a) and cat's-ears diagrams (b); s— channel and, via crossing, 
II— channel nucleon resonances (c); channel meson exchange (d). 



elements, for example, electromagnetic form factors. If 
all ingredients satisfy appropriate Ward-Takahashi iden- 
tities (WTIs), electromagnetic gauge invariance is satis- 
fied by construction. The method has been used for com- 
puting nucleon form factors both in the quark-diquark 
model [32] and from the covariant three-body equa- 
tion [33]. The scope of the approach is however more 
general: the same formalism can be exploited to derive 
generalized parton distributions from their quark-gluon 
substructure [34] , and it has been recently applied to ob- 
tain the scattering amplitudes for Nt: scattering and pion 
electroproduction at a hadronic level, i.e., from chiral ef- 
fective field theories [35, 36]. 

In the following we will use the gauging procedure to 
resolve the coupling of a hadron to two external cur- 
rents from its substructure in QCD. The hadronic scat- 
tering matrices obtained in that way are constructed 
from QCD's Green functions and hadron wave func- 
tions which, in turn, must be determined beforehand 
from their DSEs and bound-state equations. As a conse- 
quence, all hadron resonances shown in Fig. 1 can be ex- 
pressed through diagrams that involve only the degrees of 
freedom in the QCD Lagrangian. Since the method does 
not discriminate between the scattering on mesons or 
baryons, and the involved currents can have any desired 
Lorentz structure and thus describe photons as well as 
mesons, one arrives at a unified theoretical framework for 
a variety of reactions including nucleon Compton scatter- 
ing (A^7* — ?> N'^*), pion- nucleon scattering {Ntt — >■ Ntt), 
pion photo- and electroproduction {N^* — > Ntt), or tttt 
scattering. We will show in Section IV F that a rainbow- 
ladder truncation, i.e., a gluon-exchange interaction be- 
tween quark and antiquark, recovers the diagrams that 
were employed in Refs. [37, 38] to compute the tttt scat- 
tering amplitude. Our goal is to generalize that approach 
to accommodate arbitrary interaction kernels and, in par- 
ticular, to describe scattering on baryons as well. 

The manuscript is organized as follows. In Sec. II 
we provide a brief overview of the covariant bound-state 
approach. The gauging method that provides the link 
between a hadron's current and the underlying quark- 
antiquark vertices is discussed in Sec. Ill and exempli- 
fied for the derivation of hadron form factors. In Sec. IV 
we derive the relation for scattering amplitudes, exam- 



ine their ingredients in various kinematical limits and 
provide specific examples for scattering on baryons and 
mesons. We discuss possible applications and their prac- 
tical feasibility in Sec. V and conclude in Sec. VI. 



II. BOUND-STATE EQUATIONS 

In order to collect the necessary tools for the deriva- 
tion of the scattering amplitudes, we recall some ba- 
sic relations of the covariant bound-state approach in 
this section. More detailed overviews can be found in 
Refs. [24, 26, 27, 39]. The following considerations are 
general and apply for any bound state of n valence quarks 
and/or antiquarks. In practice, of course, we are primar- 
ily interested in mesons with n = 2 and baryons with 
n = 3. 

The quantity in QCD which describes a hadron micro- 
scopically is the n— quark (i.e., 2n— point) Green func- 
tion G, as well as its n— quark connected and amputated 
counterpart, the scattering matrix T that is defined via 



G — Go + Go T Go . 



(1) 



Here, Go denotes the disconnected product of n dressed 
quark propagators, e.g.: Go = S S S in the case of a 
baryon. To keep the discussion as transparent as possi- 
ble we will use a symbolic notation throughout the paper. 
All Dirac-Lorentz, color and flavor indices as well as mo- 
mentum dependencies are suppressed, and the products 
in Eq. (1) and subsequent equations are understood as 
four-momentum integrations over all internal loop mo- 
menta. 

The Green function G can be expanded in its Dyson 
series, thereby defining the n— quark kernel K that is the 
sum of m— quark irreducible components (2 < m < n): 



G = Go + Go K Go + Go K Go K Go 



(2) 



Its nonperturbative resummation yields Dyson's equa- 
tion: 



G = Go + Go G o 



G- 



Gn 



K . 



(3) 
(4) 
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FIG. 2. (Color online) Left panel: Dyson equation (6) for the scattering matrix in the three-quark system. Right panel: 
bound-state equation (9) for a baryon. 



The equivalent series for the T— matrix reads 

T = K + KGoK + ..., (5) 

with its resummed version, illustrated in Fig. 2: 

T = K + KGqT ^ (6) 
= K-^ - Go . (7) 

So far we have not gained much: we have related one 
unknown quantity (G, or T) with another {K) . The cen- 
tral observation is that hadrons must appear as poles in 
the respective n— quark Green function or, equivalently, 
in the n— quark scattering matrix. These poles will be 
generated selfconsistently upon solving Eq. (6), and they 
may be real and timelike or shifted into the timelike com- 
plex plane via a non-zero decay width. A pole in the 
scattering matrix defines a 'bound state' on its mass shell 
— —A'P. The scattering matrix at the pole assumes 
the form 

T > ^ + regular terms , (8) 

where ^ defines the hadron's bound-state amplitude and 
\1/ is its charge conjugate. The same relation holds for 
the Green function G if the amplitude is replaced by the 
covariant wave function Gq^', i.e., with quark propaga- 
tor legs attached. At the pole, upon inserting Eq. (8) 
in Dyson's equation and comparing the residues of the 
singular terms, Eq. (6) reduces to a homogeneous bound- 
state equation for the amplitude 5': 

XGo«' = * ^ T-^^ = 0, (9) 

where we have used Eq. (7) to arrive at the second ver- 
sion. Eq. (9) is the homogeneous Bethe-Salpeter equation 
for a bound state of n valence quarks and shown in Fig. 2 
for the case of a baryon. It can be solved once the kernel 
K and the dressed quark propagator S that enters Go 
have been determined. 

The total hadron momentum enters Eq. (9) as an 
external variable. The equation can be viewed as an 
eigenvalue problem for the quantity KGq with eigenval- 
ues A,(P2) [40, 41]: 

i^Go*. -A,*, ^ (10) 
T^i^-, = (Ar'-l)Go4'.. (11) 



The eigenvector describes a ground or excited-state 
hadron, with mass Mi , if its eigenvalue satisfies the con- 
dition Aj(P^ = —Mf) = 1. The ground state corresponds 
to the largest eigenvalue of the matrix KGq. 

The bound-state equation determines a hadron ampli- 
tude up to a normalization. A covariant normalization 
follows from the requirement that the T— matrix has unit 
residue at the pole, i.e. 

T > -jp^—^~py2 + 1°^^ singular terms , (12) 

where ' denotes the derivative d/dP^. Evaluating the 
relation 

{TT-^y = =^ T' = ~T{T-^)'T (13) 

at the hadron pole yields, in combination with Eqs. (8) 
and (12), the on-shell canonical normalization condition: 

*(T"^)'* = 1. (14) 

As a corollary of Eq. (11), the derivative of the ker- 
nel K that enters the normalization integral via Eq. (7) 
can be traded for the derivative of the eigenvalue on 
the mass shell, so that the normalization simplifies to: 
-A'* Go* = 1 [42, 43]. 

The numerical effort in solving a hadron's bound-state 
equation (9) is related to the Poincare-covariant struc- 
ture of the hadron amplitudes. While a pion amplitude 
^ apip, P) depends on four covariant basis elements, a nu- 
cleon amphtude '^aPjs{p,Q:P) includes already 64 [44]. 
Here, Greek subscripts are Dirac indices and p, q are 
relative four-momenta. In fact, the ability to solve the 
equations including the full covariant structure of the am- 
plitudes is a prerequisite for computing form factors and 
scattering amplitudes. While the structure of ground- 
state baryons such as the nucleon and the A-baryon is 
dominated by s waves, orbital angular momentum in 
terms of p waves, which are a consequence of Poincare 
covariance, plays an important role as well [33, 45-47]. 

The physics that generates these phenomena is en- 
coded in the dressed quark propagator S and the kernel 
K which represent the input of the equations, cf. Fig. 2. 
The quark propagator and the one-particle-irreducible 
counterpart of K are subject to the Dyson-Schwinger 
equations and, in principle, can be solved numerically 
within a suitable truncation. Moreover, they are re- 
lated by vector and axial-vector Ward-Takahashi iden- 
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FIG. 3. (Color online) Quark-antiquark vertex (16) and its behavior on the meson bound-state pole, Eq. (18). 



titles vifhich ensure electromagnetic current conserva- 
tion as well as Gell-Mann-Oakes-Renner and Goldberger- 
Treiman relations at the hadron level. 

These identities can also be utilized as a systematic 
construction principle for kernel ansatze [48]. In practice, 
the majority of calculations in this particular framework 
have been performed in the rainbow-ladder truncation, 
where the qq kernel is given by a dressed gluon-ladder 
exchange, with two bare structures for the quark-gluon 
vertices and an effective gluon propagator built from the 
gluon dressing function and a simple model for the vertex 
dressing; see [25, 49, 50] for a selection of results. The 
quark propagator can then be easily solved from its DSE. 
Similarly, solutions of a baryon's bound-state equation, 
which involves the three-quark kernel 

K = if(3) + (5""^ (g) X(2) + permutations) , (15) 

were obtained upon neglecting the ggg— irreducible part 
ii'(3), thereby defining the covariant Faddeev equation, 
and by using a rainbow-ladder ansatz for the qq contri- 
bution K(^2) [33, 44, 47]. Recent progress in constructing 
qq kernels beyond rainbow-ladder, without resorting to 
an order-by-order resummation, together with promising 
results in the meson sector [43, 51-55] provide hope that 
such expressions may eventually also be suitable for an 
implementation in baryon calculations. 

We wish to emphasize that the subsequent derivation 
is completely general and does not rely on any specific 
assumptions for the input S", Ki^2) ^-nd K(^^^ of the bound- 
state equations. We will treat these quantities essen- 
tially as black boxes and assume that they are known in 
advance. While the rainbow-ladder truncation will ul- 
timately prove useful for practical implementations, the 
derivation of the scattering amplitudes is valid for general 
kernels. 



III. HADRON CURRENTS 

A. Quark-antiquark vertices 

In order to relate the coupling of an external current, 
for instance, a pseudoscalar, vector or axialvector cur- 
rent, to the underlying description of the hadron as a 
composite object, one must specify how that current cou- 
ples to the hadron's constituents. At the microscopic 
level, the current couples to quarks. That is described 



by the quark-antiquark vertex F'*, cf. Fig. 3: 
F^ :=Go'GF^'— F^ + TGoF^, 



(16) 



where in the second step we have exploited Eq. (1). The 
quantities G, Go, T and K are those introduced in Sec- 
tion II in the quark-antiquark case. The vertex is the 
contraction of the quark-antiquark Green function with 
a 'bare' term Fq, where ^ is not necessarily a Lorentz 
index but rather a label for the type of current that will 
be identified with the 'gauging' index in Section IIIB. 
For instance, Fg can represent 75, 7^*, or 757^, with ap- 
propriate renormalization constants attached. 

Inserting Dyson's equation for G in (16) yields the in- 
homogeneous Bethe-Salpeter equation for the vertex: 



F^ = F^ 



KGnT". 



(17) 



which features the same ingredients as a meson's bound- 
state equation and thus can be solved consistently within 
the same truncation. 

Since the qq scattering matrix T contains meson 
bound-state poles, with the pole behavior from Eq. (8), 
any pole will also appear in the vertex according to (16), 
so that on the mass shell of the total momentum the ver- 
tex becomes proportional to the bound-state amplitude: 



p2__ 71 r2 



P2 + M2 



Go Tq 



P2 +M2 



(18) 



This relation holds as long as the respective meson wave 
function 5* Gq has non- vanishing overlap with the struc- 



ture Fn and the residue 



Gq Fq is thus non-zero. 



Among other consequences, this feature is the underly- 
ing reason for 'vector-meson dominance' in hadron elec- 
tromagnetic form factors: at the microscopic level, the 
coupling to a photon is represented by the quark-photon 
vertex which is a Lorentz vector. The bare term 7^^ has 
an overlap with the p— meson wave function, and that 
overlap is proportional to its electroweak decay constant. 
Consequently, the p— meson pole appears in the (trans- 
verse part of the) quark-photon vertex and transpires to 
the level of hadron form factors. If the quark-antiquark 
kernel K that enters Eqs. (9) and (17) as an input allows 
for a p — > TTTT decay, the p— meson acquires a width and 
its pole will be shifted into the complex plane. 

Eq. (18) further implies that a hadron's coupling 
to a meson can be described by the respective quark- 
antiquark vertex in the same manner as the interaction 
with a photon. Upon removing the onshell meson pole 
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together with its residue from the vertex, the couphng 
to the quark is represented by the meson's bound-state 
amphtude. This is the central feature that enables a com- 
mon description of hadron-meson and hadron-photon in- 
teractions: with a single expression for the qq kernel A', 
different types of quark-antiquark vertices can be readily 
solved from the same inhomogeneous BSE (17) which dif- 
fers only by the inhomogeneous driving term Fg . In view 
of investigating the hadron's coupling to a photon, one 
solves the BSE with the inhomogeneous term 7'^ to ob- 
tain the quark-photon vertex; for describing its coupling 
to a pion, one would implement 75 and remove the pion 
pole together with its residue from the resulting pseu- 
doscalar vertex, or solve the pion's bound-state equa- 
tion (9) directly. Apart from that, there is no conceptual 
difference in the description. Thus, whenever a quark- 
antiquark vertex appears in the following considerations, 
we will leave its Lorentz type unspecified. 



B. Gauging of equations 

While, at least in principle, the tower of QCD's Green 
functions can be solved from Dyson-Schwinger equations, 
the question of how to combine these Green functions 
for obtaining structure functions at the hadron level is 
already much less straightforward to answer. A system- 
atic construction principle to derive the coupling of a 
hadron to an external current is the 'gauging of equa- 
tions' method [29-31]. 

Diagrammatically, gauging corresponds to the cou- 
pling of an external current with given quantum numbers 
to an n— point Green function. When acting upon that 
Green function it yields an [n + 1)— point function with 
an additional leg. Gauging is formally indicated by an 
index ^ and has the properties of a derivative, i.e. it is 
linear and satisfies Leibniz' rule. At the quark level, it 
amounts to replacing each inverse dressed quark propa- 
gator by the respective quark-antiquark vertex, i.e., the 
vertex is the gauged inverse quark propagator: 



s- 



r'^^ (i + rGo)r^. (i9) 



Correspondingly, each quark propagator is replaced by 
S"'' ~ — Sr^'S", where the specific type of gauging is de- 
termined by the structure Fq. 

To derive a quark-level description of a hadron's cou- 
pling to an external current (or several external currents), 
one applies the gauging operator to the n— quark scat- 
tering matrix T. In the same way as the pole residues 
in T correspond to hadron bound-state amplitudes, the 
residue of defines a hadron's current that, depending 
on the type of Fq, includes the various form factors of 
a hadron: for example, its electromagnetic (7^), pseu- 
doscalar (75) or axial form factors (757'^). If T is gauged 
twice, the residue in T'"^ is the desired scattering ampli- 
tude that describes, for instance, Compton scattering or 
the scattering of a nucleon and a pion. With the help 



of Dyson's equation (6) and the Leibniz rule, the gaug- 
ing operation can be systematically traced back to the 
gauging of quarks and that of the kernel K. 

Before applying the formalism to determine scattering 
amplitudes, we will revisit the derivation of hadron form 
factors in the next section. We will make frequent use 
of the derivative property of the gauging operator which 
implies 



and analogous relations for other Green functions. 



(20) 



C. Derivation of form factors 

The form factors of a hadron are the Lorentz-invariant 
dressing functions of the hadron's current matrix element 
J^. The latter is a three-point function that describes 
the hadron's interaction with the external current. At 
the level of QCD's Green functions, that interaction is 
contained in the gauged n— quark scattering matrix, i.e., 
in the (2n+ 1)— point function T**. To isolate a hadron's 
form factor inside T** , one must inspect its pole structure 
that is given in terms of a pole for the incoming and one 
pole for the outgoing hadron. The pole residue defines 
the current matrix element: 



(P2 



J'^ 



(21) 



where (P/ — PiY = is the four- momentum transfer 
and, for instance in the baryon case, '^i — "^{pijqi, Pi) 
and \I'/ = "^{pf,qf,Pf) are in- and outgoing baryon am- 
plitudes with different momentum dependencies. On the 
other hand, Eqs. (20) and (8) yield 



rpfi — _ y (T^^') T 



Pf=Pf = -M^ 



(22) 



P2 + AP 



P? + A P ' 



and by comparing these two equations one obtains the 
current as the gauged inverse scattering matrix element 
between the onshell hadron amplitudes: 



(23) 



In the limit — > 0, Eq. (23) reproduces the canonical 
normalization condition: once the bound-state amplitude 
is normalized via Eq. (14), the normalization of the re- 
sulting form factors is fixed as well. If the vector WTI is 
respected throughout every stage, electromagnetic cur- 
rent conservation and thus charge conservation, e.g., for 
a pion or a proton, readily follow [32]. Similarly, if the 
axialvector WTI is correctly implemented at the level 
of the kernels, the axial form factors satisfy Goldberger- 
Treiman relations [56]. 
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To make Eq. (23) more useful for practical numerical 
implementations, we need to express the inverse scatter- 
ing matrix by the kernel K and the propagator product 
Go- Gauging Eq. (7) yields 



Go Go - K-^K''K~^ , 



(24) 



where we have denoted the gauged inverse propagator 
product by (Gq ^)'^ =: 
system it is given by 

= (s-i(g)S-i)^ = r^(8)S' 

and thus 



r''. For instance, in a two-body 
'-3-^(8)^, (25) 



Go Go = s^''s®s + s®sr^'s. 



(26) 



When acting on the onshell amplitudes, the inverse 
kernels in (24) can be replaced by Go via the bound- 
state equations K^^'^i = Go*i and ^fK~^ = ^ f Go- 
This yields the final result for the current: 



J'' = Go (r^ - K") Go 



(27) 



This relation has been exploited in numerous meson 
form-factor studies [57-60]. In rainbow-ladder trunca- 
tion: K^^s^ — 0, and the resulting meson current is ex- 
pressed by impulse-approximation diagrams alone. In 
the baryon case, the three-body kernel of Eq. (15) with 



(3) 



yields the covariant Faddeev equation. The 



gauged kernel in rainbow-ladder becomes 



(2) 



permutations , 



(28) 



which was recently implemented to compute the nu- 
cleon's electromagnetic form factors [33]. The analogue 
of Eq. (27) in the quark-diquark model was utilized 
in various quark-diquark studies of baryon form fac- 
tors [26, 32, 61, 62]. 



IV. HADRON SCATTERING 

A. Derivation of the scattering amplitude 

We now turn to the derivation of four-point functions 
within the present formalism. Generalizing Eq. (21), one 
can identify the hadron's coupling to two external cur- 
rents as the residue of the scattering matrix T^'^ that is 
gauged twice: 



(29) 



J'"' is the desired scattering matrix on the mass shell 
of the incoming and outgoing hadron. The strategy to 
obtain its quark-level decomposition proceeds along the 
same lines as for the form factors in the previous sec- 
tion. After expressing J'^'^ as a function of T, and 



T'"^, we will relate these quantities to the dressed quark 
propagator S and the kernel K as well as their gauged 
analogues. 

In the first step we use the derivative property (20) of 
the gauging operation to rewrite T'^''; 



TiT-^Y^T 



-TIT 



T 



T{T-yT 



(30) 



+ T{T-YT{T-^yT 

+ t{t-^Yt{t-^Yt 



^ T 



where the curly brackets denote a symmetrization in the 
indices and v. Upon replacing the T— matrices on the 
far left and right with 



Pi 



Af2 



and 



P? + A'P 



(31) 



respectively, we can identify the onshell residue J^'^: 



(r-i)^^T(T 



T 



(32) 



This equation is the equivalent of Eq. (23) for a four-point 
function. 

In the second step we want to relate the scattering 
amplitude J^" to the quantities K and Go. To obtain 
the first term in Eq. (32), we insert {T~^Y from Eq. (24) 
and use the relations 



Go T Go — G — Go , 
R-^TGo = GoTR-^ = G, 
R-^TR-^ = G + R-^ , 



(33) 



that follow from the definition of T in (1) and Dyson's 
equation (3). Note that the arrangement of the terms 
in the resummation is irrelevant: Gq R G — G R Go- 
Again, the inverse kernels can be replaced by Go when 
acting on the bound-state amplitudes. The result is 



[t-^Yt {t-^Y = Go [r^(G - Go) 

+ Rf'iG + R-^) R'' - R^'GT'' -Tf'G R" 



Go 



(34) 



For the second term in (32), we gauge (T ^Y from 
Eq. (24) once again: 



(T-i)"" ^ Go [r^^^ - r^'^ Go r"^}] Go 

+ R-^ [R^^'R-^R"^ - R>"']R-^ 



(35) 



Some of the terms cancel when summing up Eqs. (34) 
and (35), and the remainder is 



JA"^ = 5-^. Go 



{T - R)^^' G {T - RY 



Go^^ 



(36) 
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FIG. 4. (Color online) Upper panel: Diagrammatic representation of the scattering amplitude J^'^ , Eq. (39), in the case of 
a baryon. A symmetrization of the two currents with indices ^ and v leads to an additional crossed diagram. Lower panel: 
Decomposition of the quantity A'', cf. Eqs. (37), (41) and (44), that enters the scattering amplitude. For each of the first three 
graphs on the r.h.s. there are two further permutations with respect to the quark lines. Only the first two diagrams contribute 
in rainbow-ladder truncation. 



With the shorthand notation 

- K^" , A^'^ = r^'' - K^''' , (37) 

the final result for the current (27) and the scattering 
amplitude (36) can be written in the compact form 



J'' = *fGoA''GoVl/,, 



Go*, 



(38) 
(39) 



Eq. (39) is the central result of this work. We reiterate 
that the equation is completely general: we have not as- 
sumed any specific form for the type of currents, or the 
type of hadron, or the structure of the kernel that enters 
the equation. Fig. 4 illustrates the relation for the case 
of a baryon. We will examine its physical content in the 
next section. 



B. Ingredients of J^" and J'"' 

To analyze the features of the scattering amplitude J^"^ 
in Eq. (39), we need to examine its ingredients in more 
detail. It involves the quantities 



p/i ppi/ J^^^ j^fil- 



G, 



(40) 



i.e., the gauged disconnected propagator products, the 
gauged kernels, and the n— quark Green function G. We 
have assumed knowledge of the kernel given in terms 
of a diagrammatic expression or an ansatz, and conse- 
quently of the quark propagator S, the quark-antiquark 



vertex F'^, and the bound- state amplitude ^P. The propa- 
gator can be solved from its DSE which involves a quark- 
gluon vertex that is compatible with the expression for 
K; the vertex is solved from its inhomogeneous BSE (17); 
and the bound-state amplitudes are determined from 
their covariant bound-state equations. Thus, the remain- 
ing task is to express the quantities in Eq. (40) in terms 
of S, and K. 

We have already stated the decomposition of for 
the quark-antiquark case in Eq. (25); the generalization 
to the three-quark system is straightforward: 



F^ F^ ® S- 



2 permutations . 



(41) 



Let us now have a second look at Eq. (39). Due to the 
decomposition of G from Eq. (1) we note that the first 
term contains a sum of a disconnected contribution Gq 
and another part that involves the scattering matrix T. 
In combination with the decomposition of F'^ contained 
in A'^ via Eq. (37), the part with Go leads to diagrams in 
the baryon's scattering amplitude where the two currents 
couple to the same quark (handbag diagrams) or different 
quark lines (cat's-ears diagrams). On the other hand, the 
appearance of the scattering matrix T in the other part 
holds the promise for recovering baryon resonances in the 
s and u channels. We will return to this observation in 
more detail in Section IV C. 

Applying the gauging operation to the previous rela- 
tion for F^ yields the quantity T^'^ that enters the second 
term in (39). One obtains 

Tt"" =T'"' (g) S-^ + S-^ (gjVt"" +T^^' (g)^^ (42) 
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in the quark-antiquark case and 



currents: 



= r^'' «) s- 



2 permutations 



p{M (g, r"^ (g) 5""^ +2 permutations 



(43) 



for a system of three quarks. Here we have introduced the 
vertex T'^'^ that describes the couphng of a dressed quark 
to two external currents. We will analyze its features in 
detail in Section IV D and show that it is the source of 
those diagrams in J^'' that, on a hadronic level, amount 
to channel meson exchange. The contributions from 
the second line of Eq. (43) will provide further cat's-ears 
diagrams in the four-point function J''". 

The remaining quantities to be determined are the 
gauged kernels and K^'^. For a meson the respective 
construction is, at least in principle, straightforward: for 
a given diagrammatic expression of the quark-antiquark 
kernel K(^2), the gauged kernel is obtained by replacing all 
internal dressed quark propagators by —SV^S, and each 
further Green function by its gauged counterpart (as long 
as that object exists). The kernel in the baryon case, on 
the other hand, is the sum of two- and three-quark irre- 
ducible contributions given in Eq. (15). Therefore, the 
respective gauged kernel also picks up contributions pro- 
portional to r'" and F'"': 

K'' ^ K^^^ + (F^ (g) K(2) + perm.) 

+ (S-i (gX^^) +pcrm.) , (44) 

K^'' = K^^^ + (F^'' ® i^(2) + perm.) 

+ (F^^ ® + perm.) 

+ {S-^ (E) K^"^^ + -pci-m.) . (45) 

The object = T^' - K^, composed of Eqs. (41) 
and (44), is depicted in the lower panel of Fig. 4. 



Go A^'' 

- p2 



Go * * Go 

F2 + M2 



(46) 



where, depending on the order of the indices fi and v, 
— p2 is equal to the Mandelstam variable s or u. This is 
exactly the type of diagram illustrated in Fig. Ic. 

Eq. (39) additionally implements all ofFshell contri- 
butions of nucleon resonances through the three-quark 
Green function G. In practice, however, a fully self- 
consistent determination of G through Dyson's equa- 
tion (3) is not feasible due to its nature of being a six- 
point function. Owing to the rapidly increasing num- 
ber of basis elements and independent momentum vari- 
ables, present computational resources allow for a self- 
consistent computation of four-point functions at best, 
such as for example the nucleon bound-state amplitude, 
without truncating the covariant structure of the desired 
quantity. 

A potential alternative is to exploit Dyson's equation 
and derive an inhomogeneous Bethe-Salpeter equation 
for the five-point function 



■qi'^ := Go^GA^^Go *i 



(47) 



that enters the expression for J''''. Using G = Go 
Gq K G, one arrives at the equation 



^-r ^A'^Go 'f. + KGn *r 



(48) 



which determines VE'J' self-consistently once the driving 
term A'^Gq is known. The first term in the scattering 
amplitude (39) then acquires the form 



Go A^^Go* 



(49) 



and substituting the analogue of Eq. (48) for yields 



C. Hadron resonances 

As already anticipated above, the first term in the scat- 
tering amplitude of Eq. (39) exhibits a welcome feature: 
it involves the n— quark Green function G which contains 
intermediate hadron resonances. For instance, if the two 
currents with labels fi and v represent photons, J^^" de- 
scribes Compton scattering on a nucleon; if one of the 
two currents is of pseudoscalar nature, it describes pion 
photoproduction, and if both indices are pseudoscalar, 
one obtains the nucleon-pion scattering amplitude. In 
these cases, the Green function G involves all possible 
s— channel nucleon resonances and, via symmetrization 
in the indices fj, and v, also those in the u channel. Gom- 
parison with Eqs. (8) and (38) entails that at the re- 
spective pole locations = — of the resonances the 
scattering amplitude J^'' becomes the product of two 



^-J^GoXGo <^ 



(50) 



where the interactions with the currents were completely 
absorbed in the five-point functions ^'f j. 

Nevertheless, the attempt to solve Eq. (48) still 
presents a numerical challenge. As a first step in ex- 
ploring the features of Eq. (39), a separable pole ap- 
proximation in the spirit of Eq. (46), with intermedi- 
ate nucleon and A resonances, seems far more practi- 
cal. In that respect it will be advantageous to separate 
the 'pure' pole contribution encapsulated in T from the 
remaining non-resonant structure. Naturally, following 
such a course provides limited insight since one can only 
recover the effects of those resonances that were explic- 
itly inserted into the diagram in the beginning. In view 
of the availability of NN^, NNn, NA^ and NAt: form 
factors in the rainbow-ladder truncated Dyson-Schwinger 
approach [33, 56, 60, 63], however, such an attempt would 
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still be worthwhile as a first step. On a more general note, 
even once Eq. (48) is solved fully selfconsistently, the fi- 
delity of the rainbow-ladder truncation in view of nucleon 
resonances beyond the A(1232) is essentially unknown. 
Thus, a self-consistent determination of J^"^ should also 
be augmented by an effort to go beyond rainbow-ladder. 



D. 



and meson form factors 



In contrast to the kinematical structure of the scatter- 
ing amplitude in the s channel, the second contribution 
to Eq. (39) which involves the four-point function F^'^ is 
easier to access. F^'" describes the coupHng of a dressed 
quark to two external currents. To derive an expres- 
sion for that vertex, we start from the inhomogeneous 
BSE (17) for the quark-antiquark vertex F'^. Gauging 
that equation with an additional index v removes the 
constant inhomogeneous term and yields: 



(F^)"^ = (F;; + i^GoF'^)'' = {KGoT'^Y 
K^GoTi" +KG'^T''+KGori"' . (51) 



Since the sequence of gauging is irrelevant, the quantities 
in the second line must be symmetric under exchange of 
fj, and v. Eq. (51) is again an inhomogeneous Bethe- 
Salpeter equation, now for the quantity F'^'^, with an in- 
homogeneous term 



K^Go F^ + i^GfJF'^, 



(52) 



and can be solved selfconsistently once the quark prop- 
agator, F'^ and if are known. In comparison with the 
considerations of the previous subsection, this equation is 
relatively easy to handle since F'^'^ is merely a four-point 
function. 

Due to its structure of an inhomogeneous BSE, one can 
rearrange Eq. (51) to obtain a form analogous to (16): 

T'"' = Go^Gr'^'' = {l + TGo)r^o'' . (53) 
From Eq. (6) one has {l + TGo)K = T and thus 
PM^^ = {l+TGQ)K''GQT^' + TG'^T^' = 
= K^Go F^ - T Go A'' Go F^ . 



(54) 



The essence of this relation is that the object F'^'', just 
as the vertex F'^ itself, includes meson bound-state poles 
in the total quark-antiquark momentum which have their 
origin in the scattering matrix T. Consider, for example, 
the case where the index /i corresponds to a pion and to 
a photon. In that case, J^" from Eq. (39) describes pion 
(virtual) photoproduction, and F^ in Eq. (54) represents 
the pion bound-state amplitude The lowest-lying me- 
son pole in T corresponds to another pion, and at the 
respective mass shell of the internal pion pole F'^'' be- 
comes 



- 



P2 _^ ]^,.p 



Go A'^ Go ^ , 



(55) 



where —P^ denotes here the Mandelstam variable t and 
M the pion mass. Comparison with Eq. (38) implies that 
the residue at the pole is just the pion's electromagnetic 
current: 



P2 + M2 



(56) 



Thus, the pion's electromagnetic form factor, sketched 
in Fig. Id, will be recovered in the solution of the inho- 
mogeneous BSE (51) and enters the pion photoproduc- 
tion amplitude J^^^ through the term F^"^ which provides 
a natural description of an 'offshell pion'. In addition, 
all other meson poles whose bound-state amplitudes pro- 
duce non-zero form factors with a pion and a photon will 
appear in the vertex F'"' and subsequently in J'"' as well. 



E. Baryons in rainbow-ladder 

After having studied the general features of the scat- 
tering amplitude, we turn to practical applications and 
analyze its decomposition under specific assumptions for 
the involved kernels. The baryon's bound-state equa- 
tion (9) was recently solved for both nucleon and A 
baryons by neglecting the three-quark irreducible part 
if(3) that enters the kernel (15) and implementing a 
rainbow-ladder truncation for the quark-quark kernel 
K(2) [44, 47]. In our context, the rainbow-ladder trun- 
cation implies ~ ^ "which simplifies the structure of 
J^^'^ considerably. Here, K(2) consists of two bare quark- 
gluon vertices ^ 7'^ that are connected by an effective 
gluon propagator. In vacuum, a single photon or pion 
cannot couple to a gluon due to Furry's theorem, charge 
conjugation and flavor arguments; however, two photons 
or pions can. To account for that possibility, we retain 
the term K'f^^-^ in our following considerations. It would 
describe non-planar processes of the type (e) in Fig. 6 in 
the case of tttt scattering and yield analogous diagrams 
for the scattering on baryons. 

We simplify the notation by writing 



F'^ = F^ ®S-^®S- 



bc 



r^^Ka^r^®K^2), (57) 

be 



and similarly for F^^'' and -RT^". The index a labels 
the quark which couples to the current, or the spec- 
tator quark with respect to the two-body kernel, and 
{a, 6, c} is an even permutation of {1,2,3}. The rela- 
tions of Eqs. (41) and (43-45) reduce in rainbow-ladder 
truncation to 



^^ = ^F^ r'^"^-^(Fr + r^^r^>), 

a a 



(58) 



which yields for the quantities A'^ and A^"^ defined in 
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FIG. 5. (Color online) Contributions to the nucleon-photon Compton amplitude J*"^ in rainbow-ladder truncation, of. Eq. (63). 
All quark propagators are dressed, and each box includes further diagrams with permuted quark lines. In addition, the 
photon indices p and v must be symmetrized. T represents the quark- ant iquark and three-quark scattering matrices and 
K the two-quark irreducible kernel. The hatched nucleon amplitudes represent the quantities defined in (60). Diagram (c) 
containts channel meson exchange whereas s— and u— channel nucleon resonances appear in diagram (d). To obtain the pion 
photoproduction amplitude, replace one of the quark-photon vertices by an onshell pion amplitude; for nucleon-pion scattering, 
replace both. 



Eq. (37): 

a 

= ^ (rr (la - Ka) + " rj:> ~ k^,^) . 

a 



Similarly, to obtain the scattering matrix J^"^ , we sep- 
arate G into its disconnected part and the term that in- 
volves the T-matrix via Eq. (1): G = 1 + T . Hadron res- 
onances appear in T only, hence such a separation does 
not alter the conclusions of Section IV C. The scattering 
amplitude thus becomes 



We can strip the notation to its bare minimum by sup- 
pressing all occurrences oiGQ = S®S®Sas well: all 
amplitudes, vertices, kernels and T-matrices are ampu- 
tated and their connection via dressed quark propagators 
is implicit. Furthermore, we can absorb the two-quark 
quantities (1q — Ka) in the baryon amplitudes and de- 
fine: 



(60) 



/ J J (-1' LL i-L (Jj 



(62) 



The baryon current of Eq. (38) in rainbow-ladder trun- 
cation becomes in that condensed notation: 



a a 



(61) 



which we can rearrange according to the breakdown 
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shown in Fig. 5: 

a 

+ ^ ri^ r:> vi,f _ 1 ri'^ r:> 

-E*/rr*." (63) 

a 
aa' 
a 

Here, each hne from top to bottom matches one of the 
diagrams (a)-(e) in Fig. 5. The colored amphtudes in the 
figure correspond to and \1/ / whereas the hatched am- 
plitudes represent the quantities ^'^ from Eq. (60). 
The third line in the above equation is identical to dia- 
gram (c) since F^" from Eq. (54) simplifies in rainbow- 
ladder truncation to 



whereas the contribution from F'"^ reads 



(64) 



Meson poles in the t— channel appear in diagram (c) 
whereas s— and u— channel hadron resonances are con- 
tained in diagram (d). Handbag contributions appear in 
all diagrams except type (b) ; one can however show that 
diagrams (a), (c) and the pure handbag content of (d) 
can be merged into a single graph of type (c) where the 
quark-antiquark scattering matrix T is replaced by the 
Green function G. 



F. Mesons in rainbow-ladder 

The representation of the scattering amplitude J'"' in 
Eq. (39) is general and holds for baryons and mesons 
alike. In the meson case, a rainbow-ladder truncation 
again removes the gauged quark-antiquark kernel K^^ 
from Eqs. (38-39). Here the notation of the previous sec- 
tion is somewhat impractical and so we use the following 
abbreviations: 



(65) 



such that = 



and T'i, is defined analo- 



gously. Suppressing the notation of the propagators, the 
rainbow-ladder current becomes 



pp 



pp 



(66) 



which is just the impulse approximation where the cur- 
rent couples to the upper and lower quark lines. 

In the case of the scattering amplitude J'^", the first 
term in (39) stemming from the disconnected Green func- 
tion Go yields 



% (F^+F^){^(Ft + F;)^}vl,,; 



2f}'^ FJ^ 



(67) 



F 



p{.p.} 



(68) 



The final result. 



pfM pi^} 



Fj^Ff +F}^Ff 



f{^tf;> 

p{Mj.p^^} 



p{Mj.pj} 
p{Mj,p^} 



(69) 



is illustrated in Fig. 6 in the case of tttt scattering. 
The first line in Eq. (69) corresponds to the impulse- 
approximation diagrams of type (a) ; the second and third 
lines represent the diagrams (b) and (c) with a vertical 
and diagonal T-matrix insertion, respectively; and the 
terms with F^^" correspond to the horizontal i— channel 
T-matrix insertion, type (d), stemming from Eq. (64). 
Finally, we have retained the term K!^'^ as well which can 
appear if the gluon-exchange kernel is resolved into quark 
loops and lead to non-planar terms of type (e). How- 
ever, since the rainbow-ladder gluon is a model ansatz 
without explicit quark-loop content such a term does not 
contribute. 

Fig. 6 demonstrates that the impulse-approximation 
diagrams alone are inconsistent with the rainbow-ladder 
truncation. In Refs. [37, 38], the tttt scattering matrix 
was computed upon implementing the full set of diagrams 
(a-d). The resulting isospin amplitudes reproduce Wein- 
berg's theorem for the scattering lengths at threshold and 
the Adler zero in the chiral limit. Moreover, they show 
excellent agreement with chiral perturbation theory in 
general kinematics, and the a and p— meson poles in the 
respective s, t and u channels are readily recovered. 

In contrast to the baryon case, all diagrams in Fig. 6 
except type (e) can be computed fully self-consistently, 
with comparatively modest numerical effort, if the exter- 
nal particles are identical. This can be seen by replac- 
ing all occurrences of T in Fig. 6 by the Green function 
G via Eq. (1), whereby the impulse-approximation dia- 
grams pick up a minus sign. G satisfies Dyson's equa- 
tion (3) which, upon contracting with two pion bound- 
state amplitudes in the spirit of Eqs. (47-50), can be 
transformed into an inhomogeneous BSE: 



'^t'" ■.:^Go^G{^S^) = 
^^"^ = * 5 * + X Go ^'^'^ . 



(70) 



'i'^'^ is a four-point function with two spinor and two 
pseudoscalar legs and thus consists of 8 Dirac basis ele- 
ments. Since the two pious are onshell, it depends on four 
Lorentz-invariant kinematic variables. In fact, the result 
of Ref. [38] was obtained by a self-consistent solution of 
Eq. (70) and its subsequent implementation in the scat- 
tering amplitude J^'^ . We finally note that the same set 
of diagrams also appears in the photon four-point func- 
tion that enters the hadronic light-by-light contribution 
to the muon anomalous magnetic moment [64] . 
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(b) 



(7 




(d) 



(c) 



^ 

{? ^ 



(e) 



if 



For example: 



^ 




FIG. 6. (Color online) tttt scattering amplitude J**" in rainbow-ladder truncation, of. Eq. (69). A further symmetrization in the 
indices /i and v duplicates the number of diagrams. To obtain the pion Compton scattering amplitude, replace the amplitudes 
attached to /i and v by quark-photon vertices. 



V. DISCUSSION 

The result for the scattering amplitude in Eq. (39) pro- 
vides, in principle, a complete description of scattering 
processes in QCD if a hadron couples to two particles 
with qq quantum numbers. Thus, the approach can be 
applied to a variety of reactions such as nucleon Comp- 
ton scattering, nucleon-pion scattering, meson photo- 
and electroproduction, tttt scattering, or strangeness and 
charm production processes; but also to crossed-channcl 
reactions, for example exclusive pp annihilation into two 
photons that will be studied with the PANDA experiment 
at FAIR/GSI [65]. The diagrammatic decomposition of 
the scattering amplitude allows to correlate effects at the 
hadron level in various kinematical limits with the un- 
derlying ingredients in QCD. The appearance of the qq 
scattering matrix in the t channel and, in the case of 
baryons, the qqq scattering matrix in the s and u chan- 
nels induces the existence of bound-state poles at certain 
values of the invariant Mandelstam variables. This fea- 
ture can be exploited to compute the 'offshell behavior' 
of the pion electromagnetic form factor that appears in 
pion electroproduction [66]. Studying virtual Compton 



scattering, on the other hand, enables to investigate two- 
photon corrections to nucleon form factors [9] , and isolat- 
ing the handbag structure in DVCS can be used to draw 
conclusions about generalized parton distributions. 

The approach relegates the determination of the scat- 
tering amplitudes to the knowledge of QCD's Green func- 
tions which are encoded in the quark propagator and the 
two- and three-quark irreducible kernels. The aspiration 
to describe experimental cross sections thus becomes a 
matter of constructing a truncation of the DSEs that 
captures the physically relevant features. Naturally, in a 
first step one would implement the widely-used rainbow- 
ladder truncation for the computation of scattering am- 
plitudes. We have outlined the resulting expressions 
for scattering on baryons and mesons in Sections IV E 
and IV F, respectively, and they are illustrated in Figs. 5- 
6. We have already mentioned that a gluon-exchange in- 
teraction yields a good description of pseudoscalar and 
vector- meson properties in the light-quark range, e.g., in 
view of TTTT scattering [38]. It is, however, also well-suited 
for a broader range of hadron properties, for example: 
the charmonium and bottomonium region [67]; nucleon, 
A and masses [44, 47]; the nucleon's electromagnetic. 



13 



pseudoscalar and axial form factors [33, 56]; or the A 
electromagnetic and iVA7 and TVAtt transition form fac- 
tors [60, 63, 68, 69]. Judging from those results, its imple- 
mentation in the framework of scattering processes seems 
to be well justified. 

Nevertheless, there are important effects which are not 
captured by a rainbow-ladder interaction. One exam- 
ple is the opening of hadronic decay channels such as 
p — >■ TTTT and A — >■ Nt: that produce non-analyticities in 
the light pion-mass regime and are essential for reproduc- 
ing the characteristic cut structures in scattering ampli- 
tudes. A rainbow-ladder kernel does not provide such 
decay mechanisms but produces instead stable bound 
states. Interactions beyond rainbow-ladder are also vi- 
tal for the features of scalar, axial-vector and isosinglet 
mesons [28, 43, 52], meson radial excitations [27, 70] or 
heavy-light mesons [71]. This suggests analogous con- 
clusions when considering baryon excitations and open 
strangeness or charm production processes. In that re- 
spect it will be crucial to complement the computation 
of scattering amplitudes by a simultaneous development 
of kernels beyond rainbow-ladder truncation. 

Probably the most prominent example for missing con- 
tributions in rainbow-ladder are pion-cloud corrections 
in the chiral and low-momentum region. Their absence 
is visible in various hadron properties such as nucleon 
electromagnetic form factors [33]. Dressing the nucleon's 
'quark core' by final-state interactions with pseudoscalar 
mesons is among the goals pursued in the approaches 
listed in the introduction, such as for example by the 
EBAC group at JLAB [19]. In the Dyson-Schwinger 
framework, pion-cloud effects should ultimately be im- 
plemented at the level of QCD's Green functions. From 
a microscopic point of view, pion exchange corresponds 
to a gluon rcsummation. This allows to reconstruct the 
relevant topologies for generating phenomenological chi- 
ral cloud effects from the underlying dynamics in QCD. 
Pion-loop contributions will thus appear in all Green 
functions with non-vanishing quark content and lead to 
complicated gluon-interaction diagrams. A more man- 
ageable route was taken in Refs. [53, 64, 72], where the 
assumption that the quark-antiquark T-matrix in the 
low-energy region is dominated by pion exchange was 
used to derive a consistent truncation of the DSEs be- 
yond rainbow-ladder. Employing the resulting qq kernel 
in meson Bethe-Salpeter studies recovers desired pion- 
cloud induced features, e.g., a reduction of meson masses 
in the low quark-mass region [53]. The eventual imple- 
mentation of such hybrid kernels with quark, gluon and 
effective pion degrees of freedom in the baryon sector 
would mark a notable step forward. 

The truncation dependence of the Dyson-Schwinger 
approach finally provides also opportunities. The con- 
struction of Eq. (39) allows to trace experimentally ob- 
servable effects at the hadron level back to the underlying 
Green functions of QCD and can thus supply information 
and potential constraints on their properties. In that re- 
spect we emphasize that it is not only the behavior of 



the quark mass function that has ramifications at the 
hadron level. The quark propagator in Landau gauge 
is well determined from lattice QCD [73, 74], and even 
the rainbow-ladder truncated Dyson-Schwinger result for 
the propagator agrees reasonably well with those data. 
Equally important for hadron properties is the structure 
of the two-quark and three-quark kernels and the quark- 
gluon vertex upon which they depend. Examples are 
the question of the dominance of two- or three-quark in- 
teractions in the baryon's excitation spectrum, or the 
insufficiencies of a gluon-exchange interaction for the ob- 
servables mentioned before. Putting constraints from ex- 
periment on these (however gauge-dependent) quantities 
would be indeed desirable. 

Apart from the truncation dependence there are fur- 
ther practical difficulties in computing Eq. (39). One 
of them was already mentioned in Section IV C for the 
baryon case and concerns the computational effort in ac- 
cessing the s— channel structure in the scattering ampli- 
tude, which is the nucleon resonance region, as it involves 
the full three-quark scattering matrix. Another prob- 
lem is limited kinematical access to the phase space of 
such processes. This is due to the singularity structure 
of the quark propagator that inevitably emerges in the 
solution of the quark DSE as a consequence of analyt- 
icity. While not a problem in principle, present numer- 
ical implementations only sample the kinematically safe 
momentum ranges of the quark propagators that appear 
inside a momentum loop. This leads to limitations, e.g., 
when computing masses of excited hadrons, and restricts 
the kinematical access to hadron form factors and scat- 
tering amplitudes. For instance, accessing the nucleon 
resonance region above threshold in the or TVtt scat- 
tering amplitudes might require an extrapolation from 
spacelike photon or pion momenta. 

These considerations make clear that the computa- 
tion of scattering amplitudes can only happen in grad- 
ual steps. The primary emphasis would be put on un- 
derstanding the mechanisms in QCD that determine the 
behavior of the partial-wave amplitudes in certain kine- 
matical regions, or to analyze and disentangle resonant 
and non-resonant background effects in the s, t and u 
channels and study offshell effects. One can furthermore 
isolate the properties of the quark core by investigating 
the current-mass dependence of the results, in particu- 
lar in the region where hadrons become stable bound 
states and the pion cloud is suppressed. Eventually, with 
advanced numerics, more sophisticated interactions and 
improved kinematical coverage, a systematic description 
of scattering processes from the underlying dynamics in 
QCD seems certainly possible. 



VI. CONCLUSIONS AND OUTLOOK 

We have detailed a systematic construction to describe 
scattering processes of photons and mesons with hadrons 
in the Dyson-Schwinger framework of QCD. The ap- 
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proach is Poincare-covariant and covers both perturba- 
tive and nonperturbative momentum regions as well as 
the full quark mass range. Hadronic scattering ampli- 
tudes are resolved in terms of the underlying dynamics 
in QCD by coupling the external currents to all dressed 
quarks. Electromagnetic gauge invariance and crossing 
symmetry are therefore satisfied by construction. The 
amplitudes can be computed consistently if their micro- 
scopic ingredients, i.e., the dressed quark propagator and 
the irreducible two- and three-quark kernels, have been 
determined in advance. The resulting expressions allow 
to relate hadron resonances as well as background effects 
which appear in the s, t and u channels to the structure 
of the internal two- and three-quark scattering matrices. 

Potential applications cover a wide range of reactions 
such as Compton scattering, nucleon-pion scattering, me- 
son photo- and electroproduction or exclusive proton- 
antiproton annihilation processes. Such studies would 
further allow to access offshell effects and two-photon 
contributions in the extraction of form factors, and they 
can potentially also provide information on generalized 
parton distributions. 

We have discussed the practical feasibility of the ap- 
proach and found that, while scattering on mesons 
is comparatively unproblematic from a computational 
point of view, the application in the baryon sector will 



require substantial numerical effort. This is especially 
true for the nucleon resonance region, whereas <— channel 
effects are easier to access. We have exemplified the ap- 
proach in a rainbow-ladder truncation, where the two- 
quark kernel is represented by a dressed gluon exchange 
and the three-quark kernel is omitted. Its successful ap- 
plication in various studies of meson and baryon phe- 
nomenology suggests that this setup provides a suitable 
starting point to study hadron-meson and hadron-photon 
reactions as well. Ultimately, in view of implementing 
pion-cloud corrections and hadronic decay channels, the 
computation of scattering processes should be paralleled 
by an effort to go beyond rainbow-ladder. The derivation 
that we have presented is general and can accommodate 
such purposes. 
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